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^ , Abstract 

o 
o 

■ We study the pair production of scalar top quarks (stop, t\) in photon-photon 

collisions with the subsequent decay of the top squarks into 6-quarks and charginos 
t\ — ► bxf- We simulate this process by using PYTHIA6 for electron beam energy 
qq ■ 2E^ eam = y / s^ = 1000 GeV. A set of criteria for physical variables is proposed which 

. provides a good separation of stop signal events from top quark pair production 

which is the main background. These criteria allow us to reconstruct the mass of 
the top squark with two years of ILC running provided that the neutralino mass is 
known. 



I Introduction. 

The scalar top quark, the bosonic partner of the top quark, has attracted much 
attention as it is expected to be the lightest colored supersymmetric (SUSY) [1] particle. 

II and £r, the supersymmetric partners of the left-handed and right-handed top quarks, 
mix and the resultant two mass eigenstates t\ and t2, can have a large mass splitting. It 
is even possible that the lighter eigenstate tx could be lighter than the top quark itself [2] , 
[3]. 

Experimental searches for top squarks were performed at LEP and Tevatron. The 
present status of these experimental efforts is reflected in Figure 1, taken from [4], which 
shows that the mass of the stop is higher than 141 GeV for a mass difference between 
the stop and the lightest neutralino of about 50-70 GeV (see Fig.l). These searches 
will continue at LHC and ILC. At ILC it is planned to have the option of a photon 
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collider (PLC), as originally planned for TESLA [5]. This will be achieved by using 
backscattered photon beams by Compton scattering of laser photon beams with electron 
beams [9], [10], [11] (for recent review on this subject see [6]). It has been stressed that 




Figure 1: The current exclusion limits, obtained at Tevatron and LEP, based on the search in 
the t — ► cxi channel (taken from [4])- 



the polarization effects which appear in the interactions of backscattered laser photons 
[12]-[15] may provide additional possibilities for studying the properties of the produced 
particles (see also [5] and [7], [8]). In what follows we study the reaction 

7 + 7->*i + *i- (!) 
Among the possible ti-decay channels within the MSSM (see [16] for details), we focus 
on the decay ti — > bxf followed by the two-body chargino decay Xi — > XiW^, where one 
of the W's decays hadronically, W — > q^qj, and the other one leptonically, W — > /iz/ M [17] 1 . 
The final state of this signal process, shown in the left-hand plot of Fig.2, contains two 
6-jets and two (or more) jets (originating from the decay of one W boson) and a hard 
muon plus a neutrino (from the decay of the other W) and two neutralinos: 

77 - Ui -> bbxtx'i - bbW+W-xIxi - bbq^^XiXv (2) 

The main background process is top quark pair production with the subsequent decay 
t — > bW^ (for W's we use the same decay channels as in the stop case): 

77 -> it -> bbW + W~ -> bbqiqjfiu^. (3) 

The only difference between the final states of stop and top production (see right-hand 
plot of Fig.2) is that the stop pair production has two neutralinos which are undetectable. 
In the present paper we consider only top pair production as background. The analysis of 

1 The process e + e~ — > t\ti with the subsequent decay t — ► c\° was considered in [18] -[20]. 
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Figure 2: Left is the stop signal event diagram, Right is the top background diagram. 



the processes (2) and (3) is done on the basis of Monte Carlo samples of the corresponding 
events. Two programs PYTHIA6.1 [21] and CIRCE2 [22] were used. To simulate stop 
pair production process (1), we used the PYTHIA6.1 event generator in which the formula 
for the cross section of the stop pair production in e + e~ annihilation was replaced by the 
formula for two scalar particles (s) production 77 — > ss from [23], [14], [24], [15], (see [25] 
for the NLO corrections and [26] for more plots on differential cross sections), which takes 
into account various photon polarization states. The top background was also simulated 
with PYTHIA6.1. The program CIRCE2 was used to generate the momentum spectra of 
the backscattered photons involved in process (1). The energy of the electron beams was 
chosen to be E^ eam = 500 GeV (i.e. the total e~e~ energy is -E*- e - = v /s^T = 1000 GeV). 

In Section 2 we give the set of MSSM parameters used in our study. 

In Section 3 the important backscattered photon beam characteristics, namely, mo- 
mentum spectra and luminosity, are considered for the case of polarized photon production 
in Compton scattering of polarized laser photons and polarized electrons. 

In Section 4 we discuss some general characteristics of the signal process 77 — > t\ti 
and the main background 77 — > ti. Subsection 4.1 deals with the reconstruction of the 
invariant mass of the two-quark system stemming from the W boson decay as well as 
with the reconstruction of the invariant mass of the corresponding two-jet system. Three 
cuts, applied to select signal stop and background top events, are presented in subsection 
4.2. Further subsections include kinematic distributions for the produced stop quarks, 
6-jets, missing energy, total visible energy, and the scalar sum of the transverse momenta 
of all visible particles in the event. We always compare them in detail with those of 
top pair production. In subsection 4.6 we present the values of the signal-to-background 
ratios (S/B) which may be reached with the three cuts. We also consider the differences 
between the distributions of the kinematic variables presented in this Section for stop pair 
production and top pair production. They can be used for further separation of the signal 
events from the background if more data are available. 

In Section 5 we discuss how to discriminate between the signal muons produced in W 
boson decays and those stemming from hadron decays in the same events. 

The following two Sections deal with the invariant mass distributions. In Section 6 we 
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demonstrate that the invariant masses of all final-state hadronic jets and of all final jets 
plus the signal muon are good tools for separating the signal from the top background. 

In Section 7 we show that the invariant mass of one 6-jet and the other two non — 6-jets 
(from W decay) allows one to reconstruct the mass of the scalar top quark provided that 
the neutralino mass is known. Section 8 contains some conclusions. 



2 MSSM parameters and cross section. 

The scalar top quark system is described by the mass matrix (in the II — Ir basis) [2], 
[27] 



M? M? 

ill tLR 

Mf Mf 

tRL tRR 



(4) 



with 

M? ll = Mj + (1 - ^sin 2 Q w )cos2(3M 2 z + Mf, (5) 
M? rr = Ml + ^sin 2 Q w cos2{3M 2 z + Mf, (6) 

M L = (KrY = M Mt - foot®. (7) 
The mass eigenvalues are given by 

1 



Mf - - - 

*1,2 2 



( M L + M U =f V( M L + M U + 4 I M LJ ( § ) 



with the mixing angle 

—Mf 

cosQ- t = tLR (9) 

'\M? \ 2 + (M? -M? ) 2 

1 tLR 1 v tl t LL ' 

M? -Mf 

sinQ i = tLL 11 (10) 

'\m \ 2 + (m? — m? ' 



tLR. 1 V tl t^L ' 

In what follows we shall consider only one particular choice of the MSSM parameters 
that are defined, in the notations of PYTHIA6, in the following way: 

• M 5 = 270 GeV; 

• M^ — 270 GeV; 

• A t — —500 GeV (top trilinear coupling); 

• tanf3 = 5; 

• 11= -370 GeV; 

• Mi = 80 GeV; 

• M 2 = 160 GeV. 



Note that in PYTHIA6 Mg corresponds to Mj l (left squark mass for the third genera- 
tion) and rrifj corresponds to Mj . These parameters give M^ = 167.9 GeV, M x + = 159.2 
GeV and M x o = 80.9 GeV. This value of M^ is rather close to the mass of the top quark 
M top = 170.9 ±1.8 GeV [28]. Therefore one expects a rather large contribution from the 
top background, which means that the choice of this value of the stop mass makes the 
analysis most difficult. Finding a suitable set of cuts separating stop and top events will 
be crucial. 
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3 Photon beam characteristics. 



Let us mention two main features of photon-photon collisions. The first one is that 
the monochromaticity of the backscattered photon beams is considerably increased if the 
mean helicities A e and P c of the electron beam and the laser photon beam are chosen 
such that 2A e P c ~ —1, as has been shown theoretically in [13] 2 . In this case the relative 
number of hard photons becomes nearly twice as large in the region of the photon energy 
fractions y~i = Ej / E% eam « 0.7 — 0.85, i = 1, 2, and the luminosity in collisions of these 
photons increases by a factor of 3-4. 

The other features stems from the fact that unlike the situation at an electron-positron 
collider, the energy of the beams of backscattered photons will vary from event to event. 
As already mentioned in the Introduction, we used the program CIRCE2 [22] which gives 
the energy spectra of the colliding backscattered photons, as well as the values of photon 
beam luminosities. CIRCE2 uses as input the data files that were generated for TESLA 
using the code described in [15] and [5] and the set of beam parameters described in [29] 
and [5] 3 . We used as a reasonable approximation the CIRCE2 output spectra obtained on 
the basis of the above mentioned data files that were originally generated for E l °^ e _ = 800 
GeV and scaled them (by 1000/800) to the higher beam energy 2E^ eam = £*°* e _ = 1000 
GeV. 

The photon energy spectrum obtained in this way with CIRCE2 for this total energy 
is shown in Fig. 3 4 . Two peaks are clearly seen in this figure. The left one at a low 
photon energy is caused by multiple Compton scattering and beamstrahlung photons 
[5]. The second one, according to [13], appears in the region of hard photon production 
?/i 5 2 ~ 0.83. It shows the degree of monochromaticity of the produced backscattered 
high-energy photons. 
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Figure 3: The whole photon momentum spectrum generated by CIRCE2. 



The energy spectra and the polarizations of the photons, as generated by CIRCE2, 

2 A laser beam polarization of 100% can be assumed. An electron polarization of 85% is expected at 
the ILC. 

3 The spectra obtained by CIRCE2 are in agreement [30] with the code CAIN [31]. 
4 Examples of energy, photon polarization and 77 luminosity spectra, obtained for a set of different 
values of total energy El°_^ e _, can be seen in [13]-[15], [32]. 
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were used as input for PYTHIA for the generation of stop pair production events. Due 
to the stop pair mass threshold 2M^ , only in about 0.3% of the CIRCE2 events the 77 

energy is high enouth for the generation 25 000 77 — > ti signal events by PYTHIA. The 




Figure 4: a) correlation spectra of the energy fractions yl and y2 for all events generated by 
CIRCE2 for the photon polarizations " H — " and " — h" (J=2 case), c) correlation spectrum for 
that part of the events shown in plot a) in which the stop pair was produced, b) and d) are the 

same correlation spectra as a) and c) spectra but obtained for " + +" and " " polarizations 

(J=0 case). El°J e _ = 1000 GeV. 

correlations between the energies of two colliding photons given by CIRCE2 are shown in 
Fig. 4. The plots presented there are obtained for the choice of the backscattered photon 
beam polarization relative to the electron polarization as 2X e P c « —1, which allows an 
increase in the contribution of the colliding two-photon system having the total angular 
momentum J = (the sign can easily be changed by changing the signs of the electron 
and laser beam polarizations). 

The two-dimensional plot a) in Fig.4 shows the correlation between yi and y~2 for 
the case of all photons have J = 2 i.e., we choose the CIRCE2 parameters such that 
two colliding backscattered photons have opposite helicities 5 . Following [13], [5] and 
[26], we denote this case as " H — " and " — h". Plot b) is for J = i.e., it is for the 
choice of the CIRCE2 parameters that corresponds to the case where the two colliding 
backscattered photons have the same helicities 3 . This case is denoted as " + +" and 

5 See Fig. 2 of [32] as an illustration. 
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" ". The distribution in plot b) shows maxima at yxp ~ 0.83, which corresponds to 

the high-energy peak in Fig.3. 

The lower two plots c) and d) of Fig. 4 are 3-dimensional plots with their projections 

onto the yi y 2 plane. They also show the correlations between the energy fractions 

yi and y 2 of the backscattered photons. In these plots we include only those events that 
lead to the production of a stop-antistop it pair. The left side of Fig. 4, as well as of all 
other figures below (except Figs. 6, 7 and 22), shows the plots for the " H — " and " — h" 

polarization case (i.e., J = 2) and the right side for the " + +" and " " polarization 

case (i.e., J = 0). Plots b) and d) show the enhancement of the J = state contribution 
at 2/i j2 ~ 0.83. 

It is worth mentioning that in a real photon-photon collision experiment none of these 
cases would appear in a pure form because of the unavoidable presence of some admixture 
of other photon polarization states 6 . 
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Figure 5: The same as in Fig. 4 but with the J = 2 state contribution enhanced. 



The simultaneous change of the signs of the laser photon and beam electron helicities 
at only one side of the colliding beams keeps the sign of the product 2X e P c [12], [13] , but 
leads to a different beam configuration, which may influence the shape of the luminosity 
spectrum. 

In Fig. 5 we present the correlation plots that are analogous to those of Fig. 4 but this 
time they are for the case of the above mentioned simultaneous sign reversal of the laser 
photon and electron beam polarizations at one side of the colliding beams. It is seen from 

6 Partly this is due to the fact that the source electron beams are not 100% polarized. 
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the plots a) and c) that this combination gives an increase of the contribution of the 
two-photon system of total angular momentum J = 2. 

Finally mention the values of total photon-photon luminosities and the corresponding 
values of stop pair production cross sections (for the chosen value of the stop mass) 
obtained from CIRCE2 and PYTHIA6 for £*°' e _ = 1000 GeV: 

• for the plots shown in Fig. 4 (i.e., enhanced J = state) 

& _+ = 9.35 * 10 2 fb^year' 1 ; a hh = 2.91 fb; 
& __ = 1.02 * 10 3 fb^year- 1 ; cr' lfl = 4.96 fb. 

• for the plots shown in Fig. 5 (i.e., enhanced J = 2 state). 

LP™_ & _ + = 1.02 * 10 3 fb^year' 1 ; a* 1 * 1 = 6.13 fb; 
L™. & = 9.35 * 10 2 fb^year- 1 ; a* 1 * 1 = 4.70 fb. 



4 Distributions of kinematic variables in stop and top 
production. 

In this Section we present various plots of distributions for different physical variables 
based on 2.5 ■ 10 4 stop pair production events generated by CIRCE2 and PYTHIA6.1. 
They were weighted by the photon-photon luminosity calculated with the help of CIRCE2 
and given above for the corresponding polarizations. Analogous plots are also given for 
1.0 • 10 5 generated background top events. 

It is worth mentioning that the generation of all events, i.e. for stop and top cases, was 
done separately for both possible combinations of photon polarizations, i.e. for " + +" 

and " " and for " H " and " h" cases. In what follows we present only those 

plots which correspond to the case where the relative alignment of laser photon and beam 
electron helicities enhances the contribution of the colliding two-photon system with the 
total angular momentum J = (i.e., corresponding to Fig.4). 7 

To find the jets we used the subroutine PYCLUS of PYTHIA with the distance mea- 
sure used in the " Durham algorithm" . The parameters of this jetfinder are chosen such 
that the number of jets is exactly four. Technically, 6-jets are defined as jets that contain 
at least one B-hadron. Their decay may be identified by the presence of a secondary 
vertex [35]. 

Let us also mention that separate stop and top event samples were generated without 
and with cuts. This shows the effects of the cuts. 



4.1 Jet distributions from W decay. 

In this subsection we present some plots that were obtained without cuts. In addition 
to the two 6-jets there are also, according to the decay chain (2), two jets due to the decay 
of one W boson into two quarks W — > qi + q~j (see Fig. 2). 

Figure 6 shows the spectrum of the invariant mass Mw = M inv (quarkl + quark2) 
reconstructed from the vectorial sum of 4-momenta of these two quarks. The main features 

7 The case of J = 2 is easier for background suppression due to spin 1/2 of the top quark. 
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of these plots practically do not differ for " H — " and " — h" and the " + +" and " " 

polarization cases. We use W* to denote a virtual W boson. Plot a) is for stop pair 
production, plot b) is for top production. In plot a) of Fig. 6 one clearly sees the virtual 
nature of the W boson in the stop pair production case. Hence, in the stop case the 
invariant mass of two quarks produced in the decay of the virtual W boson (W*) is 
smaller than the mass of a real W boson. In top production (see plot b) of Fig. 6) there 
is a peak in the invariant mass distribution at the mass value of the real W boson. 
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Figure 6: The invariant mass of two quarks Myy = Mi nv (quarkl + quark2), reconstructed from 
the vectorial sum of J^-momenta of two quarks that are produced inW — > qi + q~j decay, a) stop 
pair production; b) top pair production. 




Figure 7: Number of generated events versus the reconstructed invariant mass of "all-non-b- 
jets". a) stop pair production; b) top pair production. 



Figure 7 shows the corresponding plots at the jet level. The invariant mass in the plot 
a) is built of "all-non-fe-jets" (or, shortly, "JETSw*")- One can see from plot a) that in 
the stop case the peak position is shifted to the left and a long tail for higher invariant 
masses appears. As seen from plot b), in the top case at the jet level the position of 
the W-peak remains at the same value (with a high precision) of Myy as in plot b) of 
Fig. 6, except shifting of the mean value. From comparison of plots a) and b) of Fig. 7 we 
conclude that the cut Mi nv (all — non — b — jets) < 70 GeV may allow elimination of this 
tail and a big amount of the top background. 
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4.2 The cuts. 



Starting with this subsection we shall present the signal and background distributions 
which are obtained with the following set of cuts: 
• there must be at least two 6-jets in an event: 

N b - jets > 2; (11) 

• the invariant mass of "all-non-6-jets" in an event must be smaller than 70 GeV: 

M inv (all -non-b- jets) < 70 GeV; (12) 

• the invariant mass of each 6-jet must be smaller than 10 GeV: 



M b _ jet < 10 GeV. (13) 

We note that the 6-jets produced in the top decay are much more energetic than 6-jets 
produced in the stop decay, see subsection 4.4. 



4.3 Distributions in stop events. 

In this and following subsections we present only those signal stop and background 
top events which have passed the imposed cuts as defined above. Figures 8-11 show some 
general kinematic distributions characteristic of the produced stop pair system, i.e., the 
distributions of the energy of the stop or antistop Ej^{= E stop ), its transverse momentum 
P7~ (= PT stop ), its polar angle 6^ {= THETA stop ) (all in e~e~ c.m.s.) and the invariant 

mass of the produced stop pair M inv (ti + ti)(= M inv (STOP + STOPb ar ))- In these plots 
we do not distinguish between stop and antistop. By comparing the left-hand side of these 
figures with the right-hand side one sees the effects of the different chosen polarizations 
(and corresponding luminosities). 

In Fig. 8 one can see that the stop energy spectra start close to the value of the stop 
mass Mj = 167.9 GeV. In the case of " H — " and " — h" polarizations (plot a)) the 
spectrum has a peak at E stop pa 320 GeV and it is characterized by a high mean value 
E™or> n = 311 GeV. It means that the stops produced are rather energetic. In the case of 

" + +" and " " polarizations (plot b)) the energy spectrum is softer, having the main 

peak at E stop pa 180 — 200 GeV and the mean value about 272 GeV. So, one may expect 

that the stops produced in the " ++" and " " case are on the average less energetic than 

in the " H — " and " — h" case. One can also see a second smaller peak at Estop ~ 400 
GeV. This is due to the effect of the luminosity and cross section enhancement in the 
J = case at y\ ~ z/2 ~ 0.83 (see the right-hand plots of Fig. 4). 

Figure 9 shows analogous distributions for the stop transverse momentum P7~(= 

PTgtop)- The P7~ spectrum for " + +" and " " polarizations (plot b)) is much softer 

than for " H — " and " — h" polarizations (plot a)), with mean values of 108 GeV and 214 
GeV, shown, respectively, in plots b) and a). 

The polar angle 0^ {= THETA stop ) distributions are shown in Fig. 10. One can see 
that the distribution for " H — " and " — h" polarizations (plot a)) is very different from 
that for " + +" and " " polarizations (plot b)). 
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Figure 8: Stop energy E^{= E stop ) spectra, a) " + — " and " — +" polarizations, b) " + +" and 
" " polarizations. 




Figure 9: Stop transverse momentum PJ~ (= PT s t op ) spectra, a) " H — " and " — h" polariza- 
tions, b) " + +" and " " polarizations. 




Figure 10: Stop polar angle 9(= THETA stop ) spectra, a) " H — " and " — h" polarizations, b) 
" + +" and " " polarizations. 



The invariant mass M inv (ti + ti)(= M inv (STOP + STOPbar) spectra of the produced 
stop-antistop system are shown in Fig. 11. For " H — " and " — h" polarizations (plot a)) 
it has a peak around 550 GeV, which is about 170 GeV higher than the analogous peak 
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Figure 11: Stop pair invariant mass Mi nv (t\ + ti)(= Mi nv (STOP + STOPbar)) spectra, a) 
" H — " and " — h" polarizations, b) " + +" and " " polarizations. 



at 380 GeV for " + +" and " " polarizations (plot b)). Note that the shapes of the 

distributions of the invariant mass of the stop pairs shown in Fig. 11 follow the energy 
spectra given in Fig. 8. Thus, the second peak in plot b) of Fig. 11 at M inv (ti + t\) w 800 
GeV has the same origin as the peak in the plot b) of Fig. 8 at E S top ~ 400 GeV. 

The vertical axis in the plots shows the number of stops and antistops produced per 
year (= 10 7 sec) in each bin. Taking the integral of the distributions and dividing its value 
by two (there is one stop-antistop pair in each event) one can get the total number of 
events expected per year for the applied cuts. These numbers are shown as "Integral" 
values within the statistical frames in the right upper corners of the plots. They are 
calculated by taking into account the ratio of the photon-photon luminosity in the energy 
region above the stop pair threshold over the total photon-photon luminosity. In the case 
of " H — " and " — h" polarizations this ratio is approximately 0.419. From Fig. 11 it is 

seen that the number of events per year for the " + +" and " " backscattered photon 

polarizations is equal to N ++ / = 650. It is appreciably larger than the corresponding 

number of events per year iV + _/_ + = 359 for " H — " and " — h" polarizations. 

4.4 6-quark and 6-jet distributions in stop and top production. 

In the case of stop decay into a 6-quark and a chargino, t\ — > 6xf , the jets produced 
in 6-quark hadronization are observable objects. Their features are interesting from the 
viewpoint of experimentally distinguishing the stop signal events from the top background. 

Plots a) and b) in Fig. 12 show the distributions of the energy Eb of the 6- and 6- 
quarks (we do not distinguish them in what follows) produced in the decay t\ — > bxf for 

the " H — ", " — h" and " + +", " " polarizations, respectively. Both spectra start 

at E b « 4 GeV (corresponding to the b-quark mass), and go up to E b « 37 GeV. The 
mean values are about 16 GeV and 14 GeV, respectively. Plots c) and d) of Fig. 12 are 
two analogous plots for top pair production. The mean values are E b ~ 69 GeV and 
Eb ~ 63 GeV, respectively, which is almost twice as high as the end point of the 6-jet 
energy spectra in the stop events. 

Figure 13 shows the transverse momentum PT b spectra of 6-quarks for stop (plots a) 
and b)) and top (plots c) and d)) production. Plots a) and c) are for " H — " and " — h" 
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Figure 12: b- and b-quark energy spectra, a) and b) are for the stop pair production; c) and d) 
are for the top pair production, a) and c) " H — " and " — h" polarizations, b) and d) " + +" 
and " " polarizations. 



polarizations, plots b) and d) are for " + +" and " " polarizations. Comparing the 

distributions in Figs. 12 and 13 with the corresponding ones in Figs. 8 and 9, one can see 
that in stop pair production the 6-quarks take only a small fraction of the energy and 
transverse momentum of the parent stops. 

The kinematic distributions of the 6-quarks from the top decay are quite different. The 
6-quarks produced in top decays are very energetic and have large transverse momenta. 
This difference from the stop decay is easily understandable. The stop decays into a heavy 
chargino, while the top decays into a real W boson whose mass is only half of the mass 
of the chargino M ±. Therefore, the 6-quarks produced in top decays have a larger phase 
space than the 6-quarks produced in stop decays. 

The polar angle THETAf,) distributions of the 6-quarks in stop production are 
presented in Fig. 14. Plot a) is for " H — " and " — h" polarizations, plot b) is for " ++" and 

" " polarizations. The difference between these distributions due to the polarization 

effects is clearly seen in this figure. 

In Fig. 15 the cos(b, b) distribution is shown, where cos(b, b) is the cosine of the opening 
angle between the 3- momenta of the b- and 6-quarks produced in the same stop event. Plot 

a) is for " H — " and " — h" polarizations, plot b) is for " + +" and " " polarizations. 

It demonstrates that most of the b- and 6-quarks move approximately in the opposite 
directions, but some are in the same hemisphere. Thus, in the experiment we may expect 
a similar angular distribution of the corresponding b- and b — jets in stop pair production. 
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Figure 13: b- and b-quark PT spectra, a) and b) are for the stop production; c) and d) are 

for the top case, a) and c) " H — " and " — h" polarizations, b) and d) " + +" and " " 

polarizations. 
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Figure 14: The spectra of the b- and b-quark polar angle 8b(= THETAf,) for stop production 
events, a) " H — " and " — h" polarizations, b) " + +" and " " polarizations. 



As the next step, we take into account 6-quark hadronization into a 6-jet. Figure 
16 shows the energy Ef,_j et distributions of the corresponding 6-jets. Plots a) and b) 
of Fig. 16 are for stop production and plots c) and d) are for top production. These 
and the following plots for jets were obtained using the "Durham algorithm" definition 
of the distance measure implemented in the PYCLUS jetfinder of PYTHIA. Comparing 
the upper plots a) and b) of Fig. 16 for the b- and b- jet energy distributions in stop 
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Figure 15: The spectra of cos(b,b) for stop production events, a) " +— " and" — +" polarizations, 
b) " + +" and " " polarizations. 
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Figure 16: b-jet energy spectra, a) and b) are for stop par production; c) and d) are for top 

par production, a) and c) " H — " and " — h" polarizations, b) and d) " + +" and " " 

polarizations. 



production with the upper plots a) and b) of Fig. 12 for the 6-quark and 6-quark energy 
distributions, one observes that the corresponding mean values of the jet energies in 
Fig. 16 are only about 1 GeV and 2 GeV higher than those in Fig.12. Furthermore, the 
higher end points of the energy distributions for the 6-jets and 6-jets are somewhat higher 
than for the corresponding quarks. There are not very pronounced long tails at higher 
energies and the energy spectrum begins at smaller values when passing to the jet level. 
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Figure 17: b-jet PT spectra, a) and b) are for stop pair production; c) and d) are for top 

pair production, a) and c) " H " and " h" polarizations, b) and d) " + +" and " " 

polarizations. 



Due to the different kinematics in top production, the mean values of the 6-jet and 6-jet 
energy distributions in the top case are about 6 GeV smaller than the mean values of the 
corresponding 6-quark and 6-quark energy distributions (see plots c) and d) in Fig. 12 and 
Fig.16). 

Figure 17 shows the transverse momentum PTb-j e t distributions of the 6-jets and 6-jets 
in stop production (plots a) and b)) and top production (plots c) and d)). By comparing 
with Fig. 13 for the corresponding PT distributions at the 6-quark level we can see that for 
stop production there is practically no change for " H — " and " — h" polarizations (except 
the appearance of not very large tails at high PT), and an increase in the mean values of 

about 1 GeV for " + +" and " " polarizations. In the case of the production the mean 

values of the PT distributions of the 6-jets are about 4-5 GeV smaller than those of the 
6-quarks. 

Even after application of the cut (13) the 6-jet energy and transverse momentum 
spectra in top production still remain much harder than those in the stop case. This may 
suggests to use this cut for top background reduction. Its efficiency will be demonstrated 
in subsection 4.6. 
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4.5 Distributions of the total energy deposition and scalar Pt- 
sum. 

In stop pair production the two neutralinos and the energetic neutrino from the 
W boson decay escape detection. The simulation with PYTHIA6 allows us to estimate 
the missing energy and the missing transverse momenta that are carried away by these 
particles. We also take into account the non-instrumented region around the beam pipe 
given by the polar angle intervals < 7° and O > 173°. 
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Figure 18: Missing energy E m i ss _tot distribution, a) and b) are for stop pair production; c) 
and d) are for top pair production, a) and c) " H — " and " — h" polarizations, b) and d) " + +" 
and " " polarizations. 



The distributions of the total missing energy for stop production and top production 
are presented in the upper and lower plots of Fig. 18, respectively. In stop pair production, 
see plots a) and b), the E miss _ tot spectrum starts at 190-200 GeV. In top pair production 
(plots c) and d)), where the two neutralinos are not present, the missing energy E miss ^ tot 
is much smaller. 

Figure 19 shows the distributions of the total energy deposited in the calorimeter 
("visible energy") E ca i_ tot in stop production (plots a) and b)) and in top production 
(plots c) and d)). The large missing energy in stop production (Fig. 18) is related to the 
low visible energy (Fig. 19), while in top production the low missing energy correlates with 
the large visible energy. A cut on the total visible energy of approximately E ca i^ tot < 240 
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GeV 8 would eliminate most of the top background while approximately 10% of the signal 
events are lost. 
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Figure 19: Total energy E ca \-tot distribution, a) and b) are for stop pair production; c) and 
d) are for top pair production, a) and c) " H — " and " — h" polarizations, b) and d) " + +" 
and " " polarizations. 



Another useful observable is the scalar sum of the moduli of the transverse momenta 
in an event PT sca i sum = J2i=i where the sum goes over all detectable particles in 

the event. Figure 20 shows the distributions of the scalar sum of the transverse momenta 
for stop production (plots a) and b)) and for top production (plots c) and d)). It is seen 
that the restriction PT sca i sum < 180 GeV would lead to a good separation of the stop 
signal events from the top background. 

Therefore, the two global variables E ca i_ tot and PT sca i sum could also be used in addition 
to the cuts (11)-(13), for further reduction of the top background contribution (as well as 
for reduction of the QCD background). 

4.6 Signal— to— background ratio. 

To conclude this Section, let us mention that from the "Integral" values, given in 
Figs. 11, 18-20, one may see that for the case with J = enhanced state considered here 

8 that is equivalent to setting a lower limit for the missing energy. 
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Figure 20: PT sca i sum distribution, a) and b) are for stop pair production; c) and d) are for 

top pair production, a) and c) " H — " and " — h" polarizations, b) and d) " + +" and " " 

polarizations. 



the signal-to-background ratio can be expected to be of the order S/B = 84 in the case 
of " + -" and " - +" polarizations and S/B = 140 in the case of " + +" and " - -" 
polarizations. In other words, the cuts (11)-(13) would provide a sample of selected events 
where the signal stop events dominates over the top background by about two orders of 
magnitude. The contribution of the top background is reduced by these cuts down to 
only four events per year. 

Finally, we present the efficiency values for the three cuts (11)-(13) as defined in sub- 
section 4.2. We define them as the summary efficiencies. It means that if Effi is the 
efficiency of the first cut (11), Eff 12 is the efficiency of applying the first cut (11) and then 
the second cut (12). Analogously, Effi 23 is the efficiency of the successive application of 
the cuts (11), (12) and (13). We obtained the following results, which demonstrate the 
crucial importance of the cut (13) in the reduction of the top background: 

• For SIGNAL STOP events : 

" + - +" polarizations - Eff x = 0.87 ; Eff 12 = 0.79 ; Eff 123 = 0.31 ; 

"++"&" - -" polarizations - Eff\ = 0.81 ; Eff 12 = 0.71 ; Eff 123 = 0.28 ; 
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• For BACKGROUND TOP events 



"+-"&"-+" polarizations- Effx = 0.95 ; Eff 12 = 0.13 ; S//123 = 0.00061; 
"++"&" "polarizations- = 0.95 ; E// 12 = 0.11 ; E// 123 = 0.00054. 

5 Distributions of the signal muons. 

To select the signal stop pair production events, see the left plot of Fig.2, one has to 
identify the muon from the W decay. The corresponding kinematic distributions of the 
signal muons are shown in Fig. 21 for both polarization combinations. 
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Figure 21: Energy and PT distributions of signal muons. 
polarizations, b) and d) " + +" and " " polarizations. 



a) and c) " + 



and 



There are, however, muons in the event coming from leptonic and semileptonic decays 
of hadrons. The left and right plots of Fig. 22 show the energy E^ec-mu and the transverse 
momentum PT^ ec _ mtl spectra of these muons stemming from hadron decays within the 
detector volume, for wich we took the size from TESLA TDR [5]. It can be seen that the 
decay muons have a rather small energy E^ec-mu and transverse momentum PTd ec -mu, 
their mean values are about 1.07 and 0.72 GeV, respectively. The analogous spectra 
for the signal muons in Fig. 21 show that the signal muons have a much higher energy 
E sig _ mu and transverse momentum PT sig _ mu . One can cut off most low-energy decay 
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muons rejecting those with E mu < 6 GeV. Such a cut leads to a loss of only 0.8% of the 
signal events. 
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Figure 22: Distributions muons from leptonic and semileptonic decays, a) Energy distribution; 
h)PT distribution. 



We have also studied another way to select the signal muon from W decay. If the axes 
of all four jets in the event are known, then in general the signal muon has the largest 
transverse momentum with respect to any of these jet axes. 



6 Optimal variables for further signal/background 
separation. 

Let us first consider the invariant mass of the system that contains all observable 
objects in the final state. This invariant mass is the modulus of the vectorial sum of the 
4-momenta of all four jets in an event plus the 4-momentum of the signal muon 

M inv (Alljets + n) = ^=1,2,3,4^ + P») 2 - (14) 

The distribution of this invariant mass is shown in Fig. 23. Plots a) and b) show the 
results for stop pair production while the plots c) and d) are for top pair production. As 
seen from these plots, the cut M inv (AUjets + fi) < 230 GeV will give a good separation 
of signal stop and top background events. 

An even more efficient separation of the signal and the background can be obtained 
by using the invariant mass of all jets which is the modulus of the vectorial sum of the 
4-momenta of all four jets 

M inv {Alljets) = ^(Z i= i,2,3APp 2 - ( 15 ) 

The corresponding distributions for the signal stop events (upper plots) and for the back- 
ground top events (lower plots) are shown in Fig. 24. Applying the cut M inv (AUjets) < 160 
GeV leads to a much better separation of signal stop and top background events. Both 
cuts on the invariant masses are especially useful if there is more than two year for running 
time. 
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Figure 23: Distribution of number of events versus the reconstructed invariant mass of all jets 
and signal muon Mi nv (Alljets + //). a) and b) are for stop production; c) and d) are for top 
case, a) and c) " H — " and " — h" polarizations, b) and d) " + +" and " " polarizations. 



7 Determination of the scalar top quark mass. 

We first define the invariant mass Mi nv (b-jet, JETSw)'- 9 

M inv (b — jet, JETSw) = M inv [b — jet + (all — non — b — jets)] = 

— \J (Pb~jet + Pall-non-b-jets) 2 1 (16) 

which is constructed as the modulus of the vectorial sum of the 4-momentum . . lZ . 

b-jet(b-jet) 

of the 6-jet (6-jet), plus the total 4-momentum of "all-non-6-jets" system stemming from 
the W decay (P^ u _ non _ b _ jets = Pf etlw + P^t2 w )- 

This will allow us to determine the mass of the stop quark to a good accuracy More 
precisely if the signal event contains a as the signal muon, we have to take the 

6-jet (b-jet) (see Fig.2). This is only possible if one can discriminate between b- and 6-jets 
experimentally. Methods of experimental determination of the charge of the 6-jet(6-jet) 
were developed in [34]. In this paper we use PYTHIA for the determination of b- and 
6-jets. In reality, according to [34], a 50% efficiency of the separation of 6-jets and 80% of 

9 wc follow here the notations of subsection 4.1 
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Figure 24: Distribution of number of events versus the reconstructed invariant mass of all jets 
Mi nv (All jets), a) and b) are for stop pair production; c) and d) are for top pair production. 
a) and c) " H — " and " — h" polarizations, b) and d) " + +" and " " polarizations. 



the corresponding purity can be expected. It means that to get a sample of 360 signal stop 
events one needs about two years of ILC running for collecting the necessary statistics. 

The distributions of the invariant masses of the "6-jet+(all-non-6-jets)" system in the 
case of stop pair production are shown in plots a) and b) of Fig. 25 for the two polarization 
combinations 10 . Both distributions have peaks at M inv (b-jet, JETSw*) ~ 50 GeV. It is 
worth mentioning that the mean values in both plots a) and b), as expected, do not 
depend on the polarization. 

In the case of top pair production (see the right plot of Fig. 2) the invariant mass 
distributions of the "6-jet+(all-non-(>jets)" system, given in plots c) and d) of Fig. 25, 
have peaks which approximately reproduce the input value of the top quark mass M top = 
170.9(±1.8) GeV (see Section 2). The accuracy is not high due to the strong suppression of 
the top background by the applied cuts (11)-(13). Recall that according to the "Integral" 
values of these plots, only about four top events per two years (i.e., about 1%) are expected 
to appear as an admixture to the stop signal. Taking into account this suppression of the 
top background by about two orders of magnitude, especially in the region 80 < M inv (b- 
jet, JETSw) < 140 GeV, we can conclude that the M inv (b-jet, JETSw) distributions 

10 we only used the samples of events which passed the cuts (11)-(13). 
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Figure 25: The spectra of the invariant masses Mi nv (b-jet, JETSw) of the "b-jet+ (all- non-b- 
jets)" -system, a) and b) are for stop pair production; c) and d) are for top pair production, a) 
and c) " H — " and " — h" polarizations, b) and d) " + +" and " " polarizations. 



also yield a good separation of signal and background events. 

In the stop case, one has to take into account that one of the stops decays into three 
jets plus a neutralino Xv Therefore, the right edge of the invariant mass distribution of 
the "6-jet+(all-non-6-jets)" system corresponds to the mass difference — M^o. 

A linear fit to the right slopes of the shown M inv (b-}et,JET Sw*) distributions is used 
to find their edges, i.e. the position of their intersection with the x axis. This also 
diminishes the bias due to the long tails. Plots a) and b) of Fig. 26, which show in more 
detail the right parts of plots a) and b) of Fig. 25, have a smaller binning in order to better 
define the slope regions. Our aim is to find those slope intervals which can lead to the 
reconstruction of the stop mass. We show in the plots c) and d) of Fig. 26 the chosen slope 
regions for each polarization case. The first order polynomial fit gives the intersection 
with the x-axis at Ml^ hted9e (b-}et, JETSw*) = 87 ± 0.1 GeV. Adding the mass of the 
neutralino M^o = 80.9 GeV one gets the reconstructed stop mass M~ ec ° = 167.9 ±0.1 

•M tl 

GeV which reproduces the input value Mj = 167.9 GeV. 

Note that the smallness of the error is due to the fitting of the generated Monte Carlo 
events. The value of the "Integral" given in Fig.26 shows that with two years statistics 
we may expect only about 50 selected events in the region of the slope. It would lead 
to a larger experimental error than that given above. A further Monte Carlo simulation 
taking into account detector effects would give a more realistic error. 
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Figure 26: Fitting the right-edge, a) andc):"-\ — " and" — h" polarizations, b) and d) : " ++" 

and " " polarizations, a) and b) : zoom in of the right slopes of the analogous plots in Fig. 25; 

c) and d) : the fitted regions. 



8 Conclusion. 

We have studied stop pair production in photon-photon collisions within the framework 
of the MSSM for the total energy of the e~e~ system E™ e _ = ^T e = 1000 GeV. We 
assume that the stop quark decays dominantly into a chargino and a 6-quark, ti — > bxf, 
and the chargino decays into a neutralino and a W boson, xt ~^ x\W ± , where the W 
boson is virtual. One of the two Ws decays hadronically, W + — > qq, the other one 
leptonically, W~ — > \i~v. 

The study is based on a Monte Carlo simulation with two programs. First we have 
used the program CIRCE2 which gives the luminosity and the energy spectrum of the 
colliding backscattered photon beams. The results of CIRCE2 are taken as input for 
PYTHIA6.1. This event generator is used to simulate stop pair production and decay as 
well as top pair production being the main background. 

Three cuts (11)-(13) for the signal stop events have been proposed. For the separation 
of the signal stop events from the top background events the second cut (12) and the third 
cut (13) are the most important. They restrict the value of the invariant mass of the two 
"all-non-6-jets" , produced in the decay of the virtual W boson, and the invariant mass of 
the 6-jet. This set of cuts leads to a signal-to-background ratio as large as S/B = 84 in 
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the case of the " H — " and " — h" polarizations and S/B = 140 in the case of the " + +" 

and " " polarizations. Thus, we expect about 1% admixture of top events to the stop 

signal. This is different from the more complicated situation in stop pair production at 
LHC (see, for instance, [33]). 

In addition, we have studied in detail various distributions of kinematic variables which 
can be suitable for further experimental separation of the stop signal events from the top 
background especially in the case of large statistics. For this purpose, we have worked 
out two optimal variables M inv (AUjets + /i) and M inv (AUjets), where the latter gives the 
better separation. 

We have shown that determining the end point of the distribution of the invariant 
mass M inv (b-}et,JET Sw*) of the "6-jet+(all-non-6-jets)" system allows us to measure the 
mass of the stop quark with a good accuracy based on the statistics of two years running. 
For this the mass of Xi has to be known. 

In conclusion, we can say that the 77 channel is very well suited for the study of stop 
pair production. 
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